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ABSTRACT 

We present the analysis of XMM -Newton and Swift optical-UV and X-ray observations of the Seyfert-1/QSO Mrk 509, part of an 
unprecedented multi-wavelength campaign, investigating the nuclear environment of this AGN. The XMM-Newton data are from a 
series of 10 observations of about 60 ks each, spaced from each other by about 4 days, taken in Oct-Nov 2009. During our campaign, 
Mrk 509 was also observed with Swift for a period of about 100 days, monitoring the behaviour of the source before and after the 
XMM-Newton observations. With these data we have established the continuum spectrum in the optical-UV and X-ray bands and 
investigated its variability on the timescale of our campaign with a resolution time of a few days. In order to measure and model 
the continuum as far as possible into the UV, we also made use of Hubble Space Telescope (HST) Cosmic Origin Spectrograph 
(COS) observations of Mrk 509 (part of our coordinated campaign) and of an archival Far Ultraviolet Spectroscopic Explorer (FUSE) 
observation. We have found that in addition to an X-ray power-law, the spectrum displays soft X-ray excess emission below 2 keV, 
which interestingly varies in association with the thermal optical-UV emission from the accretion disc. The change in the X-ray 
power-law component flux (albeit smaller than that of the soft excess), on the other hand, is uncorrelated to the flux variability of the 
soft X-ray excess and the disc component on the probed timescale. The results of our simultaneous broad-band spectral and timing 
analysis suggest that, on a resolution time of a few days, the soft X-ray excess of Mrk 509 is produced by the Comptonisation of 
the thermal optical-UV photons from the accretion disc by a warm (0.2 keV) optically thick (t ~ 17) corona surrounding the inner 
regions of the disc. This makes Mrk 509, with a black hole mass of about 1~3 x 10* Mq, the highest mass known system to display 
such behaviour and origin for the soft X-ray excess. 

Key words, galaxies: active - galaxies: nuclei - galaxies: Seyfert - quasars: individual: Mrk 509 - X-rays: galaxies 
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^% 1. Introduction The most prominent feature of AGN Spectral Energy 

Distributions (SEDs) is the so-called 'big-blue-bump' peaking 
in the extreme UV, which was first att r ibuted to thermal emis- 

Variability studies of Active Galactic Nuclei (AGN) over the last ^ion from the accretion disc by|Shields| ([Mi; it is thought that 

few decades have been essential in developing our understand- higher energy photons (UV) originate from small radii in the disc 

ing of the physical structure of emitting regions in the vicinity of ^here the temperature is higher, while optical photons are pro- 

the central engine. The huge amount of energy radiated away is ^uced further out. The first X-ray spectrum of a Seyfert- 1 galaxy, 

explained by the cuiTent paradigm of accretion of material onto ^^ich established the power-law shape of the X-ray continuum, 

a supermassive black hole (SMBH). The emitting regions in the ^^^ obtained for NGC 4151, using the proEortional counter on- 

vicinity of the SMBH are too small to be resolved with current ^oard the Ariel-V satellite, by| Ives et al.| (Il976|). The power-law 

telescopes and therefore clues to the processes going on in the component in the X-ray spectra of AGN is generally believed to 

centi-al regions come from analysis of their flux and spectral vari- ^e the result of Compton up-scattering of optical-UV photons 

ability. Multi-wavelength monitoring of the AGN variability is a ^om the disc by energetic electrons in a hot corona surrounding 

most eff^ective way to investigate complex phenomena such as '^- ^P^" from the power-law, the X-ray continuum often com- 

the accretion onto a SMBH and provides diagnostics of the pro- P^s^s another component, the soft X-ray excess, which was first 

cesses in the nuclear environment of AGN. discovered in the spectrum of Mrk 509 by Singh et^ (il985D 
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using data from the first mission of High Energy Astronomy 
Observatories (HEAO-1). Shortly afterwards, the presence of a 
soft exc ess was also reported in another Seyfert-1 galaxy, Mrk 
841, by Arnaud et alJ ([l985) using European X-ray Observatory 
Satellite (EXOSAT) observations. Since then, the soft excess 
has been observed in the X-ray spectra of the majority of AGN 
and its origin remains an area of active research to this date. 
Different interpretations for the nature of the soft X-ray ex- 
cess can be found in the literature such as: (1) the high en- 
ergy tail of the thermal em ission from the accretion disc (e.g. 
lArnaud et al.lll985llPounds e t al. 1986); (2) an artefact of strong, 
relativistically smeared, partially ionized absorption in a wind 
from the inner disc (e.g. Gierlihski & Done 2004) : (3) part of 
the relativisti cally blurred photoion ised d isc reflection spectrum 
modelled by iRoss & FabianI (|2005') (e.g. I Crummvet a"Lll2006l) : 
(4) 'warm' Comptonisation: up-scattering of seed disc photons 
in a Comptonising medium which has lower temperature and 
higher optical depth than the corona responsible f or the X-r ay 
power-law emission (e.g. the study by Maadziarz et al.l 19981 of 
the Seyfert-1 NGC 5548; the analyses by Middleton et al.ll2009l 
and Jin et al. 2009 of two Narrow Line Seyfert-ls (NLSls): RE 
J1034+396 and RX J0136.9-3510). 

Past multi- wavelength monitoring campaigns of Seyfert-1 
galaxies have shown striking similarities in the optical-UV and 
X-ray variability, suggesting a strong link between the emis- 
sion in these en ergy ranges, su c h as t hat provided by inverse 
Comptonisation. IWalter & FinkI ( Il993h studied the soft X-ray 
(0. 1-2.4 keV) spectra of 58 Seyfert-1 AGN, including Mrk 509, 
using /JOS'Ar observations. They reported presence of a soft X- 
ray excess above the extrapolation of the hard X-ray power-law 
in 90% of the sources. They also found that the soft X-ray ex- 
cess is well correlated to the strength of the big-blue-bump ob- 
served with the International Ultraviolet Explorer (lUE) in the 
UV band. They concluded that the big-blue-bump is an ultravio- 
let to soft X-ray feature, which has a similar shape in all Seyfert- 
1 galaxies of their sam ple. From an intensiv e multi-instrument 
campaign of NGC 4151 . lEdelson et"an (Il996h found that optical- 
UV and 1-2 keV X-ray fluxes varied together on timescale of 
days, with the X-ray flux varying with much larger amplitude 
than the optical-UV; the phase difference between UV and 1-2 
keV X-ray was consistent with zer o lag, with an upper limit of 
< 0.3 days. [Marshall et al.l (Il997h measured the spectrum and 
lightcurve of NGC 5548 over a period of 2 months using the 
Extreme Ultraviolet Explorer (EUVE) when the galaxy was also 
monitored with HST. They found the optical-UV and EUV vari- 
ations to be simultaneous with the amplitude in the EUV twice 
that in the UV. Also the shape of the EUVE spectrum was con- 
sistent with a gradual decreasing of flux from the UV through 
to the soft X-ray with no emission lines detected in the EUV 
band. Furthermore, from a monitoring study of Mrk 509 using 
Rossi X-ray Timing E xplorer (RXT E) and ground-based optical 
observations, Marshall et al.l (l2008h report a strong correlation 
between variability of the optical and hard X-ray emission on 
timescale of a few years, with optical variations leading the hard 
X-rays by about 15 days. 

The work presented here is part of a large multi-wavelength 
campaign studying the nearby Seyfert-1/QSO M rk 509. This 
AGN has a cosmological redshift of 0.034397 dHuchra et alJ 
Il993h corresponding to a luminosity distance of 145 Mpc (tak- 
ing Ho = 73 km s' Mpc"', ^a = 0.73 and Q„, = 0.27). 
The details of our multi-wavelength campaign are presented in 
[Kaastra et al. (2011). As part of this campaign. Swift monitor- 
ing was carried out before and after a series of observations by 
XMM-Newton and the International Gamma-Ray Astrophysics 



Laboratory (INTEGRAL). One of the aims of this campaign is 
determining the location and structure of the multi-phase warm 
absorber outflows in Mrk 509, which is done by monitoring 
the response of the warm absorber to intrinsic source variabil- 
ity using high resolution spectra from the Reflection Grating 
Spectrometer (RGS) onboard XMM-Newton. In this work we fo- 
cus particularly on the optical-UV and X-ray source variability 
of Mrk 509 using data from the XMM-Newton Optical Monitor 
(OM) and the European Photon Imaging Camera (EPIC-pn). 
Our goal in this paper is to try and explain the intrinsic multi- 
wavelength variability of Mrk 509 and its relation to the soft 
X-ray excess in terms of physical changes in the accretion disc 
and corona. 

The structure of this paper is as follows. In Section 2 we 
describe the observations and data reduction. Sect. 3 focuses on 
the photometric and spectroscopic corrections applied to the data 
in order to establish the continuum; the spectral modelling is 
described in detail in Sects. 4 and 5; we discuss our findings in 
Sect. 6 and give concluding remarks in Sect. 7. 

2. Observations and data analysis 

Some details of the XMM-Newton's OM and EPIC-pn observa- 
tions and of the Swift's UV/Optical Telescope (UVOT) and X- 
Ray Telescope (XRT) observations are shown in Table [T] 

2.1. OM and UVOT broad-band filters data 

The OM and UVOT data from Image-mode operations were 
taken with six broad-band filters: V, B, U, UVWl, UVM2 and 
UVW2. The OM Image mode data were processed with the SAS 
v9.0 omi chain pipeline. We performed aperture photometry on 
each image in a fully interactive way using the omsource pro- 
gram. We selected the source and background regions to extract 
the count rates, and applied all the necessary corrections, i.e. for 
the point spread function (PSF) and coincidence losses, includ- 
ing time-dependent sensitivity (TDS) corrections. The OM count 
rates were extracted from a circle of 12 pixels radius (5.8") cen- 
tred on the source nucleus. The background was extracted from 
a source-free region of the same radius. Similarly, for the UVOT, 
we used the uvotsource tool to perform photometry using the 
recommended circular aperture radius of 10 pixels (5.0") and 
applying the standard instrume ntal corrections and calibration 
according to lPoole et al. * (*2008*). For the purp ose of spectral fit- 
ting with the SPEX code ( Kaastra et al J 1 1 996h version 2.01.04j, 
we used the filter count rates and the corresponding response 
matrices, which were constructed using the effective area of the 
instruments for each of the filters. To calculate the flux density 
at the filter effective wavel ength from the count rate, we use d the 
conversion factors given in lTalavera & OMCal TeanJ (l2008h and 
Poo le et all (12008') for the OM and the UVOT respectively. Note 
that these conversion factors have been obtained from observa- 
tions of standard stars with known spectral shape; therefore the 
flux of Mrk 509 needs to be corrected for the presence of strong 
AGN emission lines in the filter bandpasses. We have thus taken 
into account the contribution to each filter bandpass of all signifi- 
cant emission lines and of the small-blue-bump (see Sect. l3.3l for 
the OM grism modelling) so that the flux at the effective wave- 
length of the filter corresponds only to the intrinsic continuum 
emission. 

Next, in order to account for the difference in the calibration 
of the OM and UVOT instruments, we performed the following 



http://www. sron.nl/spex 
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Table 1. The XMM-Newton (OM, EPIC-pn) and Swift (UVOT, XRT) observations details. 









Start time (UTC) 




Total 


exposure time at each filter (ks) 


X-ray exposure (ks)^ 


Inst. 


Obs. 


ID 


yyyy-mm-dd 


hh:mm:ss 


V 


B 


U 


Wl 


M2 


W2 


Vgrism 


EPIC-pn 


XRT 


OM 


1 


0601390201 


2009-10-15 


07:16:25 


6.0 


7.5 


7.5 


7.5 


7.5 


8.1 


5.0 


40.0 


- 


OM 


2 


0601390301 


2009-10-19 


15:47:02 


7.5 


6.0 


3.0 


7.5 


8.0 


8.7 


5.0 


44.4 


- 


OM 


3 


0601390401 


2009-10-23 


06:08:12 


6.0 


7.5 


7.5 


7.5 


7.5 


8.1 


5.0 


42.4 


- 


OM 


4 


0601390501 


2009-10-29 


07:22:07 


7.5 


7.5 


7.5 


7.5 


6.0 


4.9 


5.0 


42.4 


- 


OM 


5 


0601390601 


2009-11-02 


03:12:58 


7.5 


7.5 


7.5 


7.5 


7.9 


7.8 


5.0 


43.7 


- 


OM 


6 


0601390701 


2009-11-06 


07:27:14 


7.5 


7.5 


7.5 


7.5 


7.5 


8.5 


5.0 


43.9 


- 


OM 


7 


0601390801 


2009-11-10 


09:08:51 


7.5 


7.5 


7.5 


7.5 


7.5 


5.0 


5.0 


42.3 


- 


OM 


8 


0601390901 


2009-11-14 


08:53:43 


7.5 


7.5 


6.0 


7.5 


7.5 


6.3 


5.0 


42.3 


- 


OM 


9 


0601391001 


2009-11-18 


02:35:08 


7.5 


7.5 


7.5 


7.5 


8.9 


9.5 


5.0 


45.6 


- 


OM 


10 


0601391101 


2009-11-20 


08:07:05 


7.5 


7.5 


7.5 


7.5 


7.5 


8.2 


5.0 


43.7 


- 


UVOT 


1 


00035469005 


2009-09-04 


14:10:14 


- 


0.1 


0.1 


0.3 


- 


0.3 


- 


- 


0.9 


UVOT 


2 


00035469006 


2009-09-08 


04:06:41 


0.1 


0.1 


0.1 


0.2 


0.3 


0.5 


- 


- 


1.3 


UVOT 


3 


00035469007 


2009-09-12 


18:58:17 


- 


- 


- 


- 


0.4 


- 


- 


- 


0.4 


UVOT 


4 


00035469008 


2009-09-16 


20:48:13 


- 


- 


- 


- 


1.0 


- 


- 


- 


1.0 


UVOT 


5 


00035469009 


2009-09-20 


07:53:27 


- 


- 


- 


- 


- 


- 


- 


- 


1.1 


UVOT 


6 


00035469010 


2009-09-24 


10:47:43 


- 


- 


- 


- 


1.3 


- 


- 


- 


1.2 


UVOT 


7 


00035469011 


2009-10-02 


18:55:48 


- 


- 


- 


- 


1.0 


- 


- 


- 


1.0 


UVOT 


8 


00035469012 


2009-10-06 


00:53:57 


- 


- 


- 


- 


1.2 


- 


- 


- 


1.0 


UVOT 


9 


00035469013 


2009-10-10 


22:14:30 


- 


- 


- 


- 


1.1 


- 


- 


- 


1.1 


UVOT 


10 


00035469014 


2009-10-14 


11:05:50 


- 


- 


- 


- 


1.3 


- 


- 


- 


1.4 


UVOT 


11 


00035469015 


2009-10-18 


06:05:58 


- 


- 


- 


- 


1.0 


- 


- 


- 


0.9 


UVOT 


12 


00035469016 


2009-11-20 


07:30:53 


- 


- 


- 


1.0 


- 


- 


- 


- 


1.0 


UVOT 


13 


00035469017 


2009-11-24 


23:51:24 


- 


- 


1.0 


- 


- 


- 


- 


- 


1.0 


UVOT 


14 


00035469018 


2009-11-28 


20:10:30 


- 


- 


- 


1.2 


- 


- 


- 


- 


1.2 


UVOT 


15 


00035469019 


2009-12-02 


01:16:34 


- 


- 


- 


1.3 


- 


- 


- 


- 


1.3 


UVOT 


16 


00035469020 


2009-12-06 


07:19:20 


- 


- 


- 


1.1 


- 


- 


- 


- 


1.1 


UVOT 


17 


00035469021 


2009-12-08 


02:42:06 


- 


- 


- 


- 


- 


1.1 


- 


- 


1.1 


UVOT 


18 


00035469022 


2009-12-10 


12:44:41 


0.6 


- 


- 


- 


0.5 


- 


- 


- 


1.1 


UVOT 


19 


00035469023 


2009-12-12 


20:41:00 


- 


- 


- 


- 


0.9 


- 


- 


- 


1.0 



" Cleaned X-ray exposure times are shown. 



to normalise the UVOT fluxes against the OM for each filter. As 
shown in Table [1] we have 10 OM observations with exposures 
in each of the 6 filters. We thus obtained 10 SEDs with 6 optical- 
UV data points (see Fig.lTJ for the source, after implementing all 
the corrections described in Sect. [3] We then performed a cubic 
spUne interpolation between the data points for each SED and 
thus obtained a flux versus wavelength relation for each OM ob- 
servation. During the monitoring campaign there were two oc- 
casions when the OM and UVOT observations overlapped (or 
were very close to overlap). They correspond to OM Obs. 1 and 
UVOT Obs. 10 (filter M2), and OM Obs. 10 and UVOT Obs. 
12 (filter Wl). Using the OM SEDs of Obs. 1 and Obs. 10, we 
calculated the expected OM flux at the effective wavelengths of 
the UVOT M2 and Wl filters. Then by taking the ratio of the 
OM and UVOT fluxes, normalisation factors for the M2 and Wl 
filters of OM and UVOT were obtained. In order to calculate the 
normalisation factor for the other filters we made use of UVOT 
Obs. 2 in which all 6 filters were used. By knowing the average 
shape of the SED over the duration of the 10 OM observations 
(see Fig. [U we can scale the SED according to the normalised 
M2 and Wl fluxes of UVOT Obs. 2. These two SEDs can then be 



used to give us the normalisation factors of the remaining filters 
assuming the shape of the average OM SED remains the same 
over the observations. We find that at the same wavelengths, the 
OM fluxes are generally larger than the UVOT ones by about 
10% due to calibration differences. The time averaged SED also 
provides us with a template that can be used to calculate the flux 
at any particular wavelength for those UVOT observations miss- 
ing exposures with some of the filters, by knowing only the flux 
in one of the filters. This assumes that the shape of the aver- 
age OM SED is a good representation of the source and did not 
change during the Swift observations, which are before and af- 
ter the 37-day time span of the OM observations. So finally, the 
lightcurves in Fig. |2] show how the OM and UVOT continuum 
fluxes measured at the same wavelengths (eff'ective wavelength 
of the OM filters, shown above each panel in Fig.|2|i varied over 
the combined 100 days duration of the XMM-Newton and Swift 
observations. 
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Table 2. The FUSE and HST/COS continuum values corrected 
for Galactic extinction. The FUSE observations were taken on 
15 October 2000 and the COS observations are from 10-11 
December 2009. The FUSE values have been normalised to COS 
level. 



14.8 14.9 15.0 
log 1/ (Hz) 



Fig. 1. The OM SED data points (shown as open circles) for each 
observation. The vertical error bars which represent statistical 
and systematic uncertainties at Icr level are mostly within the 
size of the data symbols. The average spline SED (black curve) 
obtained from fitting the 10 OM observations is plotted in each 
panel for reference to show the variability. 



2.2. OM optical grism data 

The images from the OM optical grism, effectively a grating- 
prism combination, were first processed with the SAS v9.0 
omgchain pipelin e (for an overview of the OM instrument see 
iMason et al.ll200ll) . All the necessary corrections, including that 
for Modulo-8 fixed pattern noise and removal of scattered light 
features, were applied to obtain undistorted and rotated grism 
images. We then used the omg source program to interactively 
identify the zero and first dispersion order spectra of our source, 
to properly define the source and background extraction regions 
and extract the calibrated spectrum from the grism images of 
each observation. The flux calibration by default is performed 
between 3000 and 6000 A for the optical grism in the SAS soft- 
ware. However, because we do detect the Hq'/16563 emission 
line right at the long wavelength end of the first order spectrum, 
we extended the flux calibration to around 7000 A in order to 
include this emission line in our spectra. We used multiple OM 
spectra of the standard stars GD153 and HZ2 to extend the cali- 
bration by normalising them to HST spectra taken from the cal- 
SPECO database. We find the contamination from the second order 
spectrum (which overlaps the first order one at the location of 





Band 


Fa 




Instrument 


(A) 


(ergcm"^ s"' A"') 


Flux error 


FUSE 


941 - 944 


2.15 X 10-'3 


10% 


FUSE 


959-961 


1.93 X 10-'3 


10% 


FUSE 


983 - 985 


1.87 X 10-'3 


10% 


FUSE 


992 - 994 


1.72 X 10-'3 


10% 


FUSE/COS 


1165-1185 


1.68 X 10-" 


5% 


COS 


1405 - 1425 


1.47 X 10-'3 


5% 


COS 


1740 - 1760 


1.07 X 10-'3 


5% 



http://www.stsci.edu/hst/observatory/cdbs/calspec.html 



the Ha in the image) to be negligible compared with the strong 
emission from the line. 



2.3. HST/COS and archival FUSE data 

In order to measure the continuum as far possible into the 
UV band (up to the Galactic Lyman limit), we used HST/COS 
and FUSE spectra of Mrk 509. The COS observations, part of 
our coordinated Mrk 509 campaign, were taken on 10 and 11 
December 2009 over 10 HST orbits. The archival FUSE spec- 
trum of Mrk 509 is from 15 October 2000. Since both COS 
and FUSE share a common wavelength band from about 1165 
to 1185 A, the archival FUSE spectrum was easily scaled to 
the flux level at the time of the COS observations, assuming 
the continuum shape in the FUSE band did not vary. The data 
reduction proce dures and the full analysis of the COS spectra 
are reported in iKriss et aP (1201 ih . In this paper for our mod- 
elling we use measurements of the UV continuum (corrected for 
Galactic extinction) from narrow spectral bands, which are free 
of emission and absorption lines (see Table |2]i. The COS data 
processing pipeline uses a sensitivity function and the FUSE 
pipeline uses effective area curves to produce flux-calibrated 
spectra. Therefore for the purpose of fitting the UV spectra and 
the optical filter data together within the spex package, we had 
to convert the UV fluxes back to counts per second by changing 
to photon units and then multiplying by the effective area inte- 
grated over the narrow bandpasses given in Table |2l The COS 
observations were made close in time to the UVOT Obs. 18, and 
since UVOT Obs. 18 is close in flux level to that of OM Obs. 
2, we included the COS and the FUSE data when modelling the 
XMM-Newton Obs. 2 data. For all other XMM-Newton observa- 
tions, we scaled the COS and FUSE data based on the flux level 
observed in the OM optical-UV filters. 

2.4. EPIC-pn and XRT spectra 

The EPIC-pn data were reduced, starting from the ODE files, us- 
ing the standard SAS v9.0 software. During the XMM-Newton 
monitoring, the EPIC-pn camera was operating in the Small- 
Window mode with the thin-filter applied. We only selected 
events with FLAG -= 8 and #XMMEA_EP attributes and screened 
for increased activity of background particles. During the cam- 
paign no strong soft proton flares have been detected. The 
cleaned EPIC-pn exposure of each XMM-Newton observation is 
about 60 ks. Nevertheless, the EPIC-pn camera, when operat- 
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Fig. 2. Intrinsic optical-UV continuum lightcurves of Mrk 509 over 100 days at the end of 2009 at the wavelengths indicated 
above each panel. Blue squares: OM; red circles: measured UVOT; green triangles: UVOT observations for which fluxes have been 
calculated using a time-averaged SED as described in Sect. 12.11 The error bars of all lightcurves represent statistical and systematic 
calibration uncertainties at Icr level. 



ing in the Small-Window mode, is only active for about 70% of 
the time. Thus the final cleaned exposure time of each EPIC- 
pn spectrum is about 40 ks (see Table [T]i. The spectra were 
extracted from a circular region of 40" radius centred on the 
source. The background was extracted from a nearby source- 
free region of the same size on the same chip. The EPIC-pn data 
showed no evidence of significant pile-up, thus single and double 
events were selected. Response matrices were generated for each 
source spectrum using the SAS tasks arf gen and rmfgen. More 
details on the data reduction and analysis of the EPIC data are 
presented in Ponti et al. (in prep.). For the XRT data reduction, 
we used the online Swift-XRT data prod ucts generator faci lity 
with the default settings (details given in lEvans et alj|20()9l) to 
produce the X-ray spectra. 

3. Intrinsic continuum emission 

In order to establish the intrinsic continuum emission of Mrk 
509, we applied photometric and spectroscopic corrections to 
our data to take into account the processes which take place in 
our line of sight towards the central engine. Sects. [TTl 13.21 and 



1331 are apphcable to the optical-UV (OM and UVOT), and Sects. 
[331 and [33 to the X-ray (EPIC-pn and XRT). In the latter case 
the presence of Galactic and AGN warm absorption is taken into 
account during the fitting process. 



3.1. Galactic interstellar dereddening 

To correct the OM and UVOT optical-UV fluxes for interstel- 
lar reddening in our Galaxy we used the reddening curve of 
Cardelli et al. (1989), including the update for near-UV given 
bv O'Donnel]|Ql9940 . The colour exce ss £{8 - V) = 0.057 mag 
is based on calculations of lSchlegel et al. ( 1998) as shown in the 
NASA/IPAC Extragalactic Database (NED). The scalar specify- 
ing the ratio of total selective extinction Ry = AvJEiB - V) was 
fixed at 3.1. 



3.2. Host galaxy stellar emission 

To correct for the host galaxy starlight contrib ution in the OM 
and UVOT fields we have used the results of iBentz et alJ (l2009l) 



M. Mehdipour et al.: Multiwavelength campaign on Mrk 509. IV 



and iKinnev et all (Il996l) . iBentz et alJ (l2009l) have determined 
the host galaxy observed flux at the rest-frame wavelength of 
5 100 A {F I 5JQ0 pj for a sample of AGN using HST images. For 
Mrk 509, Fg^i_ gioo ^ = (2.52 ± 0.23) lO^'^ erg s"' cm-^ A"' for 
an aperture of 5.0" x 7.6". In order to calculate the host galaxy 
spectrum in the optical band, we used a model spectrum and 
scaled it based on the F j jjqq ^ value. Since the OM and UVOT 
apertures are taking in only the galaxy's innermost few kpc, the 
galaxy bulge template of Kinnev et al. ( 1996) was used. Figure[3] 
shows the host bulge model spectrum in the optical band, which 
be comes gradually neg l igible going towards the UV. 

iKotilaiiien & Ward (1 19941) have also measured the nuclear 
stellar flux contribution using different broad-band filters in a 
6"diameter circular aperture for a sample of galaxies. They have 
done this by subtracting the AGN contribution using profile fit- 
ting. The nuclear stellar flux of Mrk 509 for the standard B and 
V filters are shown in Fig. [3] which are consistent with the host 
bulge model spectrum. 

3.3. Correction for emission from the Broad Line Region 
(BLR) and Narrow Line Region (NLR) 

For each of the 10 XMM-Newton observations we have a 5 ks 
OM optical grism exposure with which we can model emis- 
sion from the BLR and NLR. Figure |4] shows the mean optical 
grism spectrum of the observations and our model. The accretion 
disc blackbody model spectrum (used in Sect. l4.2l and shown for 
e.g. in Fig. |9]l has a power-law shape in the optical band. Thus, 
to model the underlying optical continuum (3000-7000 A) in 
the OM grism spectra we used a simple power-law rather than 
disc blackbody model, which requires a wider energy band to 
be fitted properly. Our main purpose in fitting the OM grism is 
modelling the emission features and therefore a simple contin- 
uum model is preferred (see Sect. 14.21 for the modelling of the 
optical-UV continuum with the OM image data). As shown in 
Fig. |4] there is excess emission below about 4000 A which is 
likely to be the long-wavelength end of the 'small-blue-bump' 
feature, a blend of Balmer continuum and Fe ii line emission. 
Note that the low-level sinusoidal pattern below 3900 A in Fig. 
|4] is an ins trumental effect (re sidual of the Modulo-8 fixed pat- 
tern noise: iMason et al]|2001h . although the overall average flux 
level is correct. The spectral shape of the small-blue-bump is 
strongly blended and thus this complex feature was modelled 
with a cubic spline with linear spacing between 3000 and 4000 
A (in the observed frame) with 10 grid points. We modelled the 
broad and narrow emission lines using Gaussian line profiles. 
The best-fit parameters of the emission features are shown in 
Table [3] Then, having obtained a best-fit model for the mean 
spectrum we applied it to the individual observations and sub- 
sequently fitted them. This enables us to look for any changes 
in the emission lines during the XMM-Newton campaign (~ 37 
days) and see how they compare with the continuum variabil- 
ity. We do not detect any statistically significant variation in the 
lines like we observe for the continuum (see Fig.|2]). This is con- 
sistent with the results of the monitoring program of Mrk 509 by 
[Carone et al. ( 1996) who found the Hj6 i4861 and He ii i4686 
emission lines to respond to continuum variations with time lags 
of about 80 and 60 days respectively, considerably longer than 
the emission line lags measured for other Seyfert galaxies and 
also longer than the duration of the XMM-Newton campaign. 
Finally, we proceeded to correct the OM and UVOT filter data 
(described in Sect. 12.1b for the emission feature contributions in 
the different filter bands. 




4500 5000 5500 

Observed Wavelength (A) 



Fig. 3. The calculated bulge stellar spectrum of Mrk 509 host 
galaxy shown in the range observed by the OM optical grism. 
The spectrum is obtained from the HST measurement at rest 
wavelength of 5100 A given in Bentz et al. (2009) and using 
the bulge model template of Kinnev et al. (1996 ) to calculate the 
flux at other wavelengths. The two data points (green crosses) 
superimposed on the spectrum are the nuclear stellar flux for the 
standard B and V filters measured by Kotilainen & Ward, (1994); 
the horizontal bars represent the filter bandpasses. 




3000 3500 4000 4500 5000 5500 6000 6500 

Observed Wavelength (A) 

Fig. 4. The observed OM optical grism spectrum of Mrk 509 (al- 
ready corrected for Galactic extinction and the AGN host galaxy 
stellar emission) obtained by averaging the OM Vgrism spectra 
from the 10 XMM-Newton observations. The data are shown in 
black and the model in red. The prominent broad and narrow 
emission features are labelled. The low-level sinusoidal pattern 
below 3900 A is a residual effect of the Modulo-8 fixed pattern 



3.4. Galactic interstellar X-ray absorption correction 

In our modelling of the X-ray spectra the effects of the Galactic 
neutral absorption were included by applying the HOT model 
(collis ional ionisation equilibrium) in spex. Assuming LoddersI 
(l2003h abundances, the Galactic H i column density in our fine of 
sight was fixed to A^h = 4.44 x 10^° cm"^ dMurphvet al.lll996h 
and the gas temperature to 0.5 eV to mimic a neutral gas. 



3.5. Warm absorber correction 

In order to perform a multi-wavelength intrinsic variability study 
one also needs to determine the X-ray continuum before ab- 
sorption by the AGN warm absorber To this end, we used the 
SPEX XABS photoionised absorption model with the parameters of 
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Table 3. Best-fit parameters of the broad and narrow emission 
features in the mean OM optical grism spectrum. The;t'J = 1.04 
(302 d.o.f.). 



EPIC-pn Obs. 1 and Obs. 5 



Emission 


Energy flux^ 


FWHM* 


Ha ^6563 


44 ±2 


4300 ± 200 


BJ3 A4S61 


17.0 ±0.5 


4400 ± 200 


Hy ^4340 


9.0 ±0.3 


4400 ± 200 


HSAAlOl 


4.2 ±0.3 


3700 ± 300 


[Oiii]/l4959 


2.2 ±0.3 


1100 ±200 


[Oiii]/I5007 


7.8 ±0.3 


1500 ± 200 


' small-blue-bump' 


44 ±5"^ 


- 



10"'^ erg s"' ctiT^. 

km s"'. 

Calculated between 2900 and 3700 A rest-wavelength. 




Observed Energy (keV) 

Fig. 5. EPIC-pn spectra of Obs. 1 (magenta) and Obs. 5 (blue) 
binned to a minimum of 100 counts per bin. The spectrum of 
Obs. 5 has a higher flux than that of Obs. 1 at low energies. 



iDetmers et"aD (l201lh derived for the stacked 600 ks RGS spec- 
trum gathered during the 2009 campaign. Note that the correc- 
tion for absorption by the warm absorber is too small to account 
for the observed X-ray variability. The warm absorption correc- 
tion is less than 6% of the X-ray photon flux in all the observa- 
tions, whereas from Obs. 1 to Obs. 5 the photon flux increases 
by about 45%. 

4. Independent modelling of the optical-UV and 
X-ray continua 

In this section we describe the models used to fit the optical-UV 
(OM) and X-ray spectra (EPIC-pn) of Mrk 509 separately. In 
Sect. |5]we perform simultaneous broad-band fits with our final 
model. Note that all the photometric and spectroscopic correc- 
tions described in the previous section, and the cosmological 
redshift have been taken into account here. All the fitted pa- 
rameter errors quoted in this paper correspond to a A^^ of 1. 
Parameters of the models given in the tables correspond to the 
source reference frame, whereas all the figures in this paper are 
displayed in the observed frame. 

4.1. The X-ray continuum and Fe Ka line 

The EPIC-pn spectra of Obs. 1 and Obs. 5 are shown together in 
Fig. |5]to display the extremes of the X-ray variability observed 
in Mrk 509 during the XMM-Newton campaign. We began our 
spectral modelling by fitting the EPIC-pn 0.3-10 keV X-ray con- 
tinuum and the Fe Ka line at ~ 6.4 keV of Obs. 1 with a sim- 
ple power-law (POW in spex) and a Gaussian line profile (GAUS 
in spex) including the absorption models described in Sects. [3741 
and 13.51 We found that a single power-law cannot fit the con- 
tinuum well between 0.3-10 keV because of the presence of a 
soft excess below about 2 keV For a single power-law fit over 
0.3-10 keV a r of 2.27 and a xl of 7.9 (1788 d.o.f.) are ob- 
tained, whereas for a fit between 2.5-10 keV we find F = 1.79 
wdxl = 1.01 (1353 d.o.f.). Figure |6] shows the 2.5-10.0 keV 
EPIC-pn fit for Obs. 1, extrapolated to lower energies, clearly 
displaying the presence of a steeper continuum below about 2 
keV. 

We then fitted the 0.3-10 keV continuum of Obs. 1 by adding 
a second power-law component to model the soft excess, which 
resulted in a best-fit with;i'^ of 1.09 (1786 d.o.f.). The photon 
indices of the two power-laws are significantly different: F - 



Table 4. Best-fit parameters of the X-ray continuum, obtained 
with a two power-law model for the EPIC-pn spectra of each 
observation as described in Sect. 14.11 





Soft power-law: 


Hard power-law: 




Obs. 


r" 


Norm* 


F" 


Norm* 


xl 1 d.o.f. 


1 


2.75 ±0.01 


3.4 + 0.1 


1.46 + 0.02 


1.4 + 0.1 


1.09/1786 


2 


2.68 ±0.01 


3.8 + 0.1 


1.41+0.02 


1.3 + 0.1 


1.07/1822 


3 


2.75 + 0.01 


4.7 + 0.1 


1.42 + 0.02 


1.4 + 0.1 


1.22/1817 


4 


2.97 + 0.01 


3.6 + 0.1 


1.47 + 0.01 


1.4 + 0.1 


1.17/1780 


5 


2.95 + 0.01 


4.7 + 0.1 


1.46 + 0.01 


1.6 + 0.1 


1.17/1804 


6 


2.77 + 0.01 


4.5 + 0.1 


1.45 + 0.02 


1.7 + 0.1 


1.18/1838 


7 


2.75 + 0.01 


4.5 + 0.1 


1.45 + 0.02 


1.5 + 0.1 


1.17/1822 


8 


2.63 + 0.01 


4.8 + 0.1 


1.40 + 0.02 


1.4 + 0.1 


1.13/1831 


9 


2.76 + 0.01 


4.1 + 0.1 


1.42 + 0.02 


1.3 + 0.1 


1.27/1821 


10 


2.71+0.01 


4.2 + 0.1 


1.44 + 0.02 


1.6 + 0.1 


1.11/1840 



" Photon index. 

* 10" photons S-' keV"' at 1 keV 



2.75 for the soft X-ray component and F = 1.46 for the hard 
one. Figure [7] shows the two power-law best-fit to the data of 
Obs. 1 . We then applied this model to the spectra of all 10 EPIC- 
pn observations. Table |4] shows the best-fit parameters for the 
two power-law model for each observation. 

The subject of this paper is the broad-band modelling of the 
Mrk 509 continuum, so we have modelled the Fe Ka line using 
a simple Gaussian profile which gives a good fit. We omit re- 
porting parameters of this line because they are not relevant to 
the main thrust of this work, and instead refer the reader to an- 
other paper reporting our campaign (Ponti et al. in prep.), which 
is completely dedicated to the study of the Fe Ka line in Mrk 
509; the data from our long campaign make it possible to study 
the profile and variability of the Fe Ka band in great detail and 
deserve a separate publication. From previou s XMM-Newton an d 
Suzaku observations of Mrk 509 reported bv lPonti et al.l (l2009l) . 
the neutral Fe Ka line does not show signs of relativistic reflec- 
tion and the ionised part of the line is complex. Furthermore, in 
another forthcoming paper (Petrucci et al. in prep.) the analy- 
sis will include the INTEGRAL data and a reflection model will 
be applied to fit the Mrk 509 spectra, as higher energy data can 
constrain the reflection parameters more tightly. 
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EPIC-pn Obs. 1 




0.5 1 2 5 10 

Observed Energy (keV) 

Fig. 6. Obs. 1 EPIC-pn single power-law fit (including a 
Gaussian Fe Ka line, Galactic and warm absorption) over the 
2.5-10 keV range. The fit is extrapolated to lower energies, dis- 
playing the presence of a soft excess below about 2 keV. The 
data are shown in black and model in red. Residuals of the fit, 
(Observed-Model)/Model, are displayed in the bottom panel. 

EPIC-pn Obs. 1 




0.5 1 2 

Observed Energy (keV) 



10 



Fig. 7. Obs. 1 EPIC-pn two power-law fit (including a Gaussian 
Fe Ka line. Galactic and warm absorption) over the 0.3-10 keV 
range. The data are shown in black and model in red. Residuals 
of the fit, (Observed-Model)/Model, are displayed in the bot- 
tom panel. 



Next, we examined the Swift XRT spectra of Mrk 509 to get 
a measure of the X-ray flux before and after the XMM-Newton 
campaign. However, one needs to note that the XRT exposures 
are very short at around 1 ks (see Table[TJ, and also the effective 
area of the XRT is very much lower than that of EPIC-pn; thus 
the XRT spectra are of very low quality, and make it impossible 
to perform any accurate spectral modelling. Therefore, we omit 
reporting parameters of XRT fits and instead in Fig. [8] we show 
how the soft X-ray (0.3 keV) and hard X-ray (4 keV) contin- 
uum fluxes varied for both XRT and EPIC-pn over the 100 days 
duration of our campaign. 



4.2. The optical-UV continuum 

Apart from line emission from the BLR and NLR (which 
has been taken into account in our modelling of the optical 
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Fig. 8. X-ray continuum flux lightcurves of Mrk 509 over 100 
days at 0.3 keV (top panel) and 4.0 keV (bottom panel). The red 
circles (with large error bars) represent the XRT points and the 
blue squares represent the EPIC-pn points. The bandwidths used 
to calculate the 0.3 keV fluxes ai-e 0.02 keV for EPIC-pn and 0.3 
keV for XRT, and the bandwidths used to calculate the 4.0 keV 
fluxes are 0.1 keV for EPIC-pn and 2.0 keV for XRT. 



continuum of Mrk 509 - see Sect. 13.3b . the dominant fea- 
ture of the optical-UV spectra of AGN is the big-blue-bump: 
this was first attributed to thermal emission from the accre- 
tion disc by Shields (1978), which has since become the gen- 
erally agreed paradigm. Annuli of the disc at different radii 
have different temperatures and therefore the disc spectrum is 
the weighted sum of different blackbody spectra. To fit the big- 
blue-bump of Mrk 509, we used the disc blackbody model 
in SPEX (DBB) which is based on a geometri cally thin, opti- 
cally thick Shakura-Sunyaev accretion disc ( Sha kura & SunvaevI 
Il973l) . The radial temperature profile of the accretion disc is 
given by T{r) ^ {3GMM[l - iRi„/ry/^]/(^nr^o-)}^''^ where T 
is the local temperature at radius r in the disc, G the grav- 
itational constant, M the mass of the SMBH, M the accre- 
tion rate, cr the Stefan-Boltzmann constant and /?;„ is the ra- 
dius at the inner edge of the disc. The flux at frequency v 
from the disc seen by an observer at distance d is given by 

Fy = [4-Khv^ COS i/(c^d^)] ff ""' (e'"'/*^^) - 1)-' rdr where / is the 
inclination angle of the disc, /?out the radius at the outer edge 
of the disc, h the Pla nck constant and c the speed of light (e.g. 
see lFrank et al.ll2002l) . The parameter /?out cannot be constrained 
by spectral fitting since most of the optical-UV radiation orig- 
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Fig. 9. The best-fit model derived in Sect. |4] independently for the optical-UV and X-ray bands for each XMM-Newton observation. 
The absorption-corrected EPIC-pn data are show in green and the black lines are the two power-law model fits to the EPIC-pn 
including a Gaussian for modelling the Fe Ka line. The red circles represent the OM data and the blue lines are the disc blackbody 
models. In panel 2, the purple squares are the HST/COS data and the brown diamonds are the scaled FUSE data. The average model 
obtained from the 10 XMM-Newton observations is shown in each panel (in grey) for reference to show the variability. 
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Table 5. Best-fit parameters of the optical-UV continuum, ob- 
tained from a disc blackbody model fit to the OM data of each 
observation as described in Sect. 14.21 



UV (2120 A) vs soft X-ray (0.3 keV) count rote 



Obs. 


J a 


Norm * 


xllA.o.i. 


1 


3.82 ± 0.06 


6.9 ± 0.6 


0.6/11 


2 


3.89 + 0.06 


6.7 ± 0.6 


0.6/11 


3 


3.94 + 0.07 


6.6 ±0.6 


0.6/11 


4 


4.08 ± 0.08 


6.2 ±0.6 


0.8/11 


5 


4.15 ±0.08 


6.0 ±0.6 


0.8/11 


6 


4.00 ± 0.07 


6.5 ±0.6 


0.6/11 


7 


3.95 ± 0.07 


6.6 ±0.6 


0.6/11 


8 


3.86 ± 0.06 


6.9 ± 0.6 


0.6/11 


9 


3.94 ± 0.07 


6.5 ±0.6 


0.7/11 


10 


3.86 ± 0.06 


6.9 ± 0.6 


0.6/11 



° Maximum temperature in the disc in eV. 

* Normalisation area (A = R- cos /) in 10^^ cm^ 

^ in ^ 



inates from the inner regions of the disc for a SMBH, thus we 
fixed /?out to 10''7?in. The fitted parameters of the model are the 
normalisation A - R^ cos / and the maximum temperature of the 
disc r^ax, which occurs at r = (49/36)/?in. 

Table |5] shows the best-fit parameters of the disc blackbody 
model for each XMM-Newton observation. The Tmax increases 
from its minimum in Obs. 1 to its maximum in Obs. 5 and then 
gradually decreases back to the level of Obs. 1 towards the end of 
the XMM-Newton campaign. The normalisation area of the disc 
blackbody remains unchanged within errors, which implies that 
7?in and inclination / parameters have not significantly varied. 
Thus the variability in the disc temperature comes from a change 
in the accretion rate. Figure|9]shows the OM SED data points and 
their corresponding disc blackbody best-fit models, and also the 
EPIC-pn two power-law best-fit models of Sect. 14.11 for all the 
10 XMM-Newton observations. 

So far we have modelled the optical-UV and X-ray data in- 
dependently of each other; in Sect. |5] we make simultaneous fits 
to both data sets in order to establish a model that accounts for 
the broad-band spectrum of Mrk 509 and for its observed vari- 
ability. 



4.3. Correlation between optlcal-UV and X-ray count rates 

Figure [8] (blue square data points in the top panel) indicates 
that the soft X-ray continuum flux of Mrk 509 varied over the 
10 XMM-Newton observations in a very similar fashion to the 
optical-UV continuum (blue square data points in Fig.|2]). In or- 
der to explore the link between the optical-UV and X-ray emis- 
sion in a model-independent way, we looked at the correlation 
between their count rates. Fig.fTOlshows the UV (2120 A) count 
rate plotted against the soft X-ray (0.3 keV) count rate for the 10 
XMM-Newton observations. In fact this correlation is as strong 
as the one between the UV and optical count rates (the Pearson 
correlation coefficient between the count rates in the UV and op- 
tical filters is about 0.85, which is statistically significant with 
probability of no correlation, from the Student's t-test, given by 
p = 0.002). Furthermore, the correlation between the optical- 
UV and X-ray count rates becomes smaller as the X-ray energy 
gets larger; this is shown in Fig.[TT] So the correlation is indeed 
strongest at those X-ray energies where the soft excess emission 
is strongest (below ~ 1 keV), implying a strong link between the 




0.9 1.0 1.1 

Soft X-ray counts/s 



1.3 



Fig. 10. The UV (2120 A) count rate from the OM UVW2 filter 
plotted versus the EPIC-pn soft X-ray (0.3 keV) count rate for 
the 10 XMM-Newton observations. The bandwidth used for the 
0.3keVbinis0.02keV. 
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Fig. 11. The correlation between the optical/UV count rates from 
the 6 OM filters and the EPIC-pn X-ray count rate for the 10 
XMM-Newton observations at energies from 0.3 keV to 10 keV. 
Black asterisks: V (5430 A) filter, black diamonds: B (4500 A), 
orange triangles: U (3440 A), filled green circles: UVWl (2910 
A), unfilled red circles: UVM2 (2310 A), blue squares: UVW2 
(2120A). The Pearson correlation coefficients are shown on the 
left y-axis and the associated null hypothesis probability p (i.e. 
probability of no correlation) from Student's t-test on the right 
y-axis. The dashed black line indicates p at 5%, which shows 
that below ~ 1 keV the correlation is statistically significant at 
95% confidence. 



optical-UV and the soft X-ray excess emission. We investigate 
this further in the next section by applying a broad-band model 
to all the data. 



5. Broad-band modelling of the continuum using 
Comptonisation 

The results of the previous section suggest that the continua 
in the optical-UV and soft X-ray bands may be physically re- 
lated, i.e. they may be the results of an underlying process link- 
ing the two. Comptonisation up-scattering the optical-UV pho- 
tons to X-ray energies could be such a process. Thus we ap- 
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Table 6. Best-fit parameters of the broad-band warm Comptonisation modelling described in Sect. |5]for each XMM-Newton obser- 
vation. 





Disc blackbody: 




Comptonisation: 




Power-law: 




Obs. 


^ max ^ 


Norm* 


Tseed^ 


Tel 


r' 


Norm'^ 


r« 


Norm^ 


;r;/d.o.f. 


1 


1.96 ±0.07 


3.9 + 0.8 


1.96(c) 


0.21 + 0.01 


16.5 + 0.4 


1.1 + 0.1 


1.89 + 0.01 


3.7 + 0.1 


1.24/1796 


2 


2.02 ± 0.07 


3.7 + 0.8 


2.02 (c) 


0.22 + 0.01 


16.2 + 0.3 


1.2 + 0.1 


1.90 + 0.01 


3.9 + 0.1 


1.19/1832 


3 


2.16 + 0.07 


3.1+0.7 


2.16(c) 


0.19 + 0.01 


17.7 + 0.3 


1.4 + 0.1 


1.96 + 0.01 


4.8 + 0.1 


1.33/1827 


4 


2.26 + 0.09 


2.8 + 0.7 


2.26 (c) 


0.19 + 0.01 


17.8 + 0.4 


1.3 + 0.1 


1.95 + 0.01 


3.7 + 0.1 


1.34/1790 


5 


2.36 + 0.09 


2.5 + 0.8 


2.36 (c) 


0.19 + 0.01 


18.2 + 0.4 


1.7 + 0.1 


1.98 + 0.01 


4.6 + 0.1 


1.33/1814 


6 


2.13 + 0.08 


3.4 + 0.8 


2.13(c) 


0.19 + 0.01 


17.7 + 0.4 


1.4 + 0.1 


1.94 + 0.01 


4.8 + 0.1 


1.31/1848 


7 


2.07 ± 0.07 


3.6 + 0.7 


2.07 (c) 


0.19 + 0.01 


17.6 + 0.4 


1.3 + 0.1 


1.95 + 0.01 


4.7 + 0.1 


1.32/1832 


8 


2.01 ± 0.07 


3.9 + 0.8 


2.01 (c) 


0.21+0.01 


16.7 + 0.3 


1.3 + 0.1 


1.93 + 0.01 


4.9 + 0.1 


1.24/1841 


9 


2.18 + 0.08 


2.9 + 0.7 


2.18(c) 


0.21+0.01 


16.5 + 0.4 


1.4 + 0.1 


1.92 + 0.01 


4.1 + 0.1 


1.34/1831 


10 


1.99 + 0.07 


3.9 + 0.8 


1.99(c) 


0.20 + 0.01 


17.1 + 0.4 


1.2 + 0.1 


1.89 + 0.01 


4.5 + 0.1 


1.27/1850 



° Maximum temperature in the disc in eV. 

* Normalisation area (A = R- cos /) in 10^' cm^. 

^ in ^ 

"-' Temperature of the seed photons in eV (coupled to T^^,^). 
'' Plasma temperature in keV. 
*■ Optical depth. 



f 10^^ photons s"' keV 
* Photon index. 



at 1 keV. 



plied the COMT component in spex to fit the optical-UV and X- 
ray data simultaneously. C OMT is based on the Comptonisation 
model of iTitarchukI d 19941) . with improved approximations for 
the parameter yS(T) which characterises the photon distribution 
over the number of scatterings which the soft photons undergo 
befo re escaping the plas ma; for more details see the spex man- 
ual (iKaastra et al.ll2010h . The seed photons in this component 
have a Wien law spectrum; we coupled its temperature Tseed to 
the disc temperature Tmax which was left as a free parameter in 
the broad-band fitting. The up-scattering Comptonising plasma 
was chosen to have a disc geometry; its parameters are the tem- 
perature Te and the optical depth r. The COMT component as- 
sumes spherical symmetry for the flux calculation. As for our 
modelling in Sect.|4] all the absorption and emission corrections 
discussed in Sect. [3] and the cosmological redshift, were imple- 
mented in our modelling here. In addition to the Comptonisation 
component for modelling the optical-UV continuum and the soft 
X-ray excess, a power-law component and a Gaussian line were 
also included to model the hard X-ray continuum and the Fe 
Ka Une; the power-law was smoothly broken at low energies 
before overshooting the optical-UV flux of the disc blackbody 
and becoming unphysical. In a forthcoming paper on our cam- 
paign (Petrucci et al. in prep.), higher energy INTEGRAL data 
will also be included in the broad-band spectral fitting, and the 
hard X-ray continuum will be modelled by Comptonisation of 
the disc photons in a hot corona, instead of using a power-law 
model as done here. In this paper, however, we focus on the di- 
rect relation between the disc blackbody emission and the soft 
X-ray excess and we prefer to use an adequate but simpler rep- 
resentation of the higher energy part of the spectrum. Table |6] 
shows the best-fit parameters of our broad-band model for each 
XMM-Newton observation. 

Figure [12] shows the broad-band Comptonisation model for 
XMM-Newton Obs. 2, displaying individual components con- 
tributing to the total model spectrum. Figure [13] shows luminos- 
ity of the X-ray power-law (top panel) and luminosity of the 
Comptonisation component (bottom panel) plotted versus the 
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Fig. 12. Best-fit broad-band model for XMM-Newton Obs. 2 as 
described in Sect. [5] The red circles represent the OM data, 
the purple squares the HST/COS data, the brown diamonds are 
the scaled FUSE data and in green are the absorption-corrected 
EPIC-pn data. The grey line is the disc blackbody component, 
the dashed magenta is the warm Comptonisation component and 
the olive colour line represents the broken power-law and in- 
cludes a Gaussian for modelling the Fe Ko- line. The black line 
is the total model spectrum. 



disc blackbody bolometric luminosity for the 10 XMM-Newton 
observations. Interestingly, there is a strong linear correlation be- 
tween the soft excess modelled by warm Comptonisation and the 
disc emission, whereas there is no such obvious correlation be- 
tween the power-law and the disc emission. 

6. Discussion 

6. 7 . Broad-band continuum variability of Mrk 509 

Optical-UV continuum lightcurves of Mrk 509 (Fig. |2]l display 
very similar variability in the 6 filters during the 100 days total 
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law is +0.09 (which is not significant with p - 0.81). Also 
the correlation coefficient between the X-ray luminosity of the 
Comptonised soft excess and the power-law is +0.17 {p - 0.63), 
which suggests that the fluxes from the two components are not 
linked over the probed timescales. But there is a significant cor- 
relation with coefficient of +0.88 {p = \0'^) between the power- 
law photon index and the disc blackbody luminosity. This be- 
haviour is similar to NGC 7469 (Nandra et al. 2000) in which 
the UV and hard X-ray fluxes were uncorrected (like in Mrk 
509), but the UV flux and the hard X-ray photon index were 
correlated. Fur t hermo re, on much larger timescales (few years), 
[Marshall et al.' (2008) have found a correlation between optical 
and RXTE hard X-ray data. In general, correlations between dif- 
ferent energy bands depend on the timescales and resolutions of 
the monitoring observations. 



1.4 



1.2 



1.0 



4.5 5.0 5.5 
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Fig. 13. Top panel: the best-fit X-ray power-law luminosity cal- 
culated between 0.2-10 keV, plotted versus the disc blackbody 
bolometric luminosity for the 10 XMM-Newton observations. 
Bottom panel: 0.2-10 keV luminosity of the soft excess mod- 
elled by warm Comptonisation, plotted versus the disc black- 
body bolometric luminosity for the 10 XMM-Newton observa- 
tions. The luminosities are obtained from the broad-band mod- 
elling described in Sect. |5] The blue line in the bottom panel 
is a fit to the luminosities with linear model: Ljoft exces,s - 
0.448 Ldisc - 1.007, where the luminosities are in units of 
erg s" 



duration of the Swift and XMM-Newton campaigns. Specifically, 
during the 10 XMM-Newton observations, the optical-UV con- 
tinuum smoothly increased from its minimum in Obs. 1 to its 
maximum in Obs. 5 (18 days after Obs. 1) and then gradu- 
ally decreased returning to almost the level of Obs. 1 towards 
the end of the XMM-Newton campaign. In Sect. 14. 31 we showed 
that there is a strong model-independent correlation between the 
optical-UV and soft X-ray flux, whereas there is a lack of cor- 
relation between the optical-UV and hard X-ray data (see Figs. 
fTOlandfTTTi. Similarly, using our broad-band model. Fig. [13] in- 
dicates a strong linear correlation between the luminosity of the 
disc blackbody and that of the Comptonised soft X-ray excess: 
^soft excess - 0.448 Ldisc " 1.007, whcrc the luminosities are 
in units of 10"^^ erg s '. However, there is no obvious correla- 
tion between the luminosity of the disc and that of the X-ray 
power-law component. The computed Pearson correlation coef- 
ficient between the luminosity of the disc and the Comptonised 
soft excess component is +0.97 (which is statistically significant 
with p = 10"^), and between the disc and the X-ray power- 



6.2. The origin of the soft X-ray excess in Mrk 509 

Here we discuss diflferent interpretations for the origin of the soft 
X-ray excess, and discuss why Comptonisation is the most likely 
explanation for its presence in Mrk 509. 

One of the early explanations for the soft X-ray excess emis- 
sion in AGN was its identification with the hig h energy tail 
of ffiermal emission fro m the accretion disc (e.g. lArnaud et al.l 
ll985llPoundsetal ." 1986). As shown in Fig|9] it is evident that 
in the case of Mrk 509 the disc blackbody emission from a ge- 
ometrically thin, optically thick Shakura-Sunyaev disc, cannot 
extend to the soft X-ray band. Often it is not possible to measure 
the peak of the disc emission for AGN as it falls in the EUV 
gap. The temperature of the disc scales with the mass of the 
black hole as M"''''*, and so more m assive black hol es tend to 
have lower disc temperatures (see e.g. lPetersonlll997h . The disc 
blackbody emission for Mrk 509, which has a relatively high 
black hole mass of 1-3 xlO'^ Mq (see Sect. 16.3b has a low disc 
temperature, which causes the peak of ffie big-blue-bump to be 
in the detection range of the HST/COS. 

Another interpretation of t h e soft X-ray excess in AGN 
is given by ' Gierlihski & Dona (12004 '). who suggest it could 
be an artefact of absorption by a relativistically smeared, par- 
tially ionised disc wind. They analysed publicly available XMM- 
Newton EPIC spectra of 26 radio-quiet Palomar-Green (PG) 
quasars and initially modelled the soft X-ray excess with a 
Comptonisation component. They found a similar temperature 
(0.1 to 0.2 keV) for the whole sample, which they claim to be 
puzzling since the objects in their sample have a large range in 
mass and luminosity, and hence disc temperature. For this reason 
they suggest the soft X-ray excess may not be real emission, and 
propose an alternative solution based on absorption. The rela- 
tivistically smeared ionised absorption was modelled with xstar 
(iBautista & Kallmanll200lh. convolved with a Gau ssian veloc- 
ity dispersion of ~ 0.2 c. lGierhnski & Pond (1200 4*) suggest that 
the smearing is such that it ehminates ffie possibility of detect- 
ing any absorption lines in the RGS spectrum. The absorption 
creates a wide dip in ffie EPIC spectrum around 1 keV which re- 
sults in an apparent soft excess at lower energies. However, we 
note that by fitting the soft X-ray excess using Comptonisation, 
there is inevitably large uncertainty in constraining the param- 
eters of the plasma when using only EPIC data as there is de- 
generacy in the fitting process. In fact, excluding the optical-UV 
data we can get a good fit to the EPIC data of Mrk 509 with a 
range of values for the plasma parameters from Tg - 0.2 keV and 
T = 18 to Te = 1.0 keV and t = 6. Therefore, broad-band mod- 
elling of simultaneous optical to X-ray data including the peak 
of the big-blue-bump emission (like performed in this work) is 
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essential to constrain the parameters of the Comptonisation. In 
most AGN, however, the peak falls in the EUV band where there 
are no data available, which causes more uncertainty when mod- 
ell ing. The small scatter in t he parameters of the plasma found 
bv lGierliriski & Donel(l2004D is likely to be due to the small en- 
ergy window used for modelling: only using broad-band obser- 
vation s can show the real scatter. Moreover, iGierliriski & Dona 
(|2004|) have strongly constrained their model for all the AGN 
in their sample by fixing the disc blackbody temperature at 10 
eV, the temperature of the hot component at 100 keV and the 
photon index of the component used to fit the soft X-ray band 
at 2. These constraints have strong impact on the best-fit results. 
Furthermore, variability of the soft X-ray excess is crucial in 
understanding its origin; from our broad-band variability study 
of Mrk 509, using both EPIC and OM data (together with the 
HST/COS and FUSE UV measurements), we have found the soft 
X-ray to be 'real' emission, which varies in a similar fashion to 
the disc emission. 

The soft X-ray excess has also been associated with the 
relativistically blurred photoionised disc reflection model of 
[Ross & Fabian (2005). According to this model, the X-ray 
power-law illuminates a relativistic accretion disc and produces 
a reflection spectrum which includes a Compton 'hump' at about 
30 keV, a strong fluorescent Fe Ka line at 6.4 keV, and a soft X- 
ray excess, composed of many emission lines that are blurred by 
relativistic motion in the accretion disc. Crummy et al. (2006) 
have reported that the reflection model makes a better fit to 
the EPIC-pn spectra of 22 Type-1 PG quasars and 12 Seyfert-1 
galaxies than the conventional model of power- law and black- 
body. However, no optical-UV data are fitted in ICrummv et al.l 
(2006) and so their reflection modelling is not tested at a broad- 
band level. Furthermore, fitting the soft X-ray excess using this 
reflection model requires fine-tuning of the ionisation parameter 
and inclination angle of the disc, and almost maxirnally ro tating 
black holes for all the objects. Done & NavakshinI (|2007|) have 
explored further the disc reflection model to check whether it 
can account for the soft X-ray excess in AGN, using the addi- 
tional constraint of hydrostatic balance on the structure of the 
illuminated disc atmosphere. They conclude that reflection from 
a hydrostatic disc cannot produce the soft X-ray excess, and con- 
stant density disc models require fine-tuning and suppression of 
the intrinsic continuum to produce the largest observed soft ex- 
cesses. 

For the case of Mrk 509, we find that the significantly vari- 
able part of the X-ray continuum is the soft excess component, 
which has strong correlation with the variability of the optical- 
UV disc emission. The hard X-ray power-law displays smaller 
variability than the soft excess and seems unconnected to the 
other components over the probed timescales. The reflection 
spectrum is produced as a result of the illumination of the disc 
by the power-law and the blurred reflection component is gener- 
ally expected to be less variable than the power-law. Thus it is 
unlikely that in Mrk 509 the observed soft X-ray excess, which 
varies in a very similar fashion to the disc intrinsic emission, is 
caused by reflection over the probed timescales. In the forthcom- 
ing paper by Petrucci et al. (in prep.), where the analysis of our 
Mrk 509 campaign will be extended to higher energy by the use 
of INTEGRAL data, blurred reflection will be tested in broad- 
band modelling of the combined XMM-Newton and INTEGRAL 
data sets. 

The more likely possibility for the origin of the soft X-ray 
excess in Mrk 509 is up-scattering of the disc photons in a warm 
Comptonising corona with lower temperature and higher optical 
depth than the one responsible for the X-ray power-law (i.e. the 



hot corona). From broad-band spectra l analysis of the Seyfert- 
1 NGC 5548, 'Magdzi arz et"an (Il998h have found that the soft 
excess requires a separate continuum component which can be 
fitted by Comptonisation of thermal photons from a cold disc 
{Tm-dx ~ 3.2 eV), in a warm (~ 0.3 keV), optically thick (r ~ 30) 
plasma. The authors also find the optical-UV and soft X-ray 
fluxes in NGC 5548 to be closely correlated. Their resu lts are 
similar to t hose found for Mrk 509 in this work. Middlet on et al.l 
(l2009h and iJin et all 1^09) also report that the soft X-ray ex- 
cess in the NLSls RE J1034H-396 and RX J0136.9-3510 can 
be explained as the result o f warm Comptonis ation of the disc 
emission. Furthermore, Done & Kubota (2006) have shown that 
an optically-thick Comptonising corona over the inner regions of 
the disc in some Black Hole Binaries (BHBs, e.g. XTE J1550- 
564) can distort a standard disc spectrum to produce a strong, 
st eep tail extending to h igher energies (shown in Figs. 4 and 5 
of iDone & Kubotall2006 ): this is seen in the Very High State of 
BHBs, when the X-ray spectra show both a strong disc compo- 
nent and a strong high-energy tail (although in this state of BHBs 
the hard X-ray power-law seems to be usually absent). The high- 
energy tail in the Very High State of BHBs may be analogous in 
nature (although the temperature of the Comptonised gas is sub- 
stantially higher) to the tail of the warm Comptonisation which 
appears as a so ft excess in the X-ray spectra of Seyferts like 
Mrk 509. Also, IZdziarski et al.1 (1200 ll) have fitted the spectrum 
of the BHB GRS 1915-H105 in the Gamma-Ray State, by warm 
Comptonisation of the disc blackbody in a plasma with a tem- 
perature of 3.6 keV and an optical depth of 4.4 (see their Fig. 
3b). 

Our broad-band variability study of Mrk 509 has shown the 
existence of a strong link between the disc emission and the soft 
X-ray excess, suggesting that the disc photons are up-scattered 
in a warm Comptonising corona close to the inner disc, with a 
temperature of about 0.2 keV and an optical depth of about 17 
(see Table|6]l. The scattered fraction of the original disc photons 
by the corona is C/[l - exp(-T)], where C/ is the covering frac- 
tion of the corona over the disc and t is the optical depth of the 
corona; this estimation neglects any angle dependencies of the 
seed and scattered photons. For Mrk 509, C/ of the warm corona 
is estimated to be about 0.25. Furthermore, we calculate that the 
fraction of accretion power dissipated into the warm corona is 
also about 0.25. 



6.3. Black hole mass and accretion rate in Mrk 509 

From the disc blackbody fits we can estimate the accretion rate 
M for Mrk 509 using the normalisation A = R^ cos / and the tem- 
perature r, = \3GM M/(?,7tR^^o-)V/'^ To do this the mass of the 
black hole is required. lPeterson et alj (12004 ) have calculated the 
black hole mass M of 35 AGN from previously published broad 
line reverberation-mapping data: M - /cTcentO"fine/G, where 
/ = 5.5 is the scaling factor, Tcent is the emission line time lag 
relative to the continuum, characterised by the centroid of the 
cross correlation function, crjine is the width of the emission line, 
c the speed of light and G the gravitational constant. For Mrk 509 
the black hole mass is reported to be 1.43 x 10** Mq. However, 
here we re-calculate the black hole mass using Tcent as measured 
bv Carone et al. ( 1996) : this is a paper entirely dedicated to the 
determination of Tcent for the case of Mrk 509 where Tcent is 
found to be 80.2 days for the HyS line. We also use the H/3 line 
width measured in this work: crii„e = FWHM/ Vln 256 - 1869 
km s"' . Adopting these values we calculate a black hole mass of 
3.0 X 10^ Mq for Mrk 509, which is twice as large as the value 



14 



M. Mehdipour et al.: Multiwavelength campaign on Mrk 509. IV 



reported in lPeterson et alJ (l2004l) . So for a Schwarzschild geom- 
etry black hole, from the normalisation A and temperature T,, 
we find M has a range of 0.24 to 0.34 Mq y"' during our XMM- 
Newton campaign. Assuming a reasonable inclination angle / of 
30° to 45°, we obtain a range of 6.0 to 7.0 Ra for the inner-disc 
radius R,,,. 



7. Conclusions 

In order to study the intrinsic optical-UV and X-ray continuum 
of the Seyfert-1/QSO hybrid Mrk 509 (using XMM-Newton OM 
and EPlC-pn, Swift UVOT and XRT, HST/COS and FUSE), we 
have taken into account various processes occurring in our line 
of sight towards the central source, such as Galactic reddening 
and absorption, host galaxy emission, BLR and NLR emission 
and the AGN multi-phase warm absorption. We conclude that: 

1. The optical-UV continuum of Mrk 509 is consistent with 
thermal emission from a geometrically thin, optically thick 
cold accretion disc with a maximum disc temperature of 
about 2 eV. The X-ray continuum can be represented by two 
distinct components: a power-law with F ~ 1.9 and a soft 
X-ray excess below 2 keV 

2. We have investigated the variability of the optical-UV and 
X-ray continuum over the 100 days duration of our monitor- 
ing campaign with a resolution of a few days. We have found 
that the variability of the soft X-ray excess is very similar in 
profile to that of the optical-UV emission from the disc. The 
flux variability of the higher energy X-ray power-law com- 
ponent (albeit smaller), on the other hand, is uncorrected 
with the variability of the soft X-ray excess and that of the 
disc emission over the probed timescales. There is however a 
correlation between the power-law photon index and the disc 
blackbody flux. 

3. We can account for the origin of the soft X-ray excess and its 
variability by Comptonisation of the disc photons in a warm 
(about 0.2 keV), optically thick (r ~ 17) corona surrounding 
the inner disc. The X-ray power-law is likely to originate 
from Comptonisation in a separate optically-thin hot corona, 
possibly located further out, in the form of an atmosphere 
surrounding the disc. 

4. Explaining the soft X-ray excess using warm 
Comptonisation has previously been suggested for the 
Seyfert-1 NGC 5548 and the NLSls: RE J1034-I-396 and 
RX J0136. 9-3510. This interpretation is analogous to the 
inner-disc warm-corona model which explains the existence 
of the high-energy tail observed in some BHBs spectra in 
the Very High State. Mrk 509 is the highest mass black 
hole system known to display such variability and origin 
for the soft excess. Only variability studies based on multi- 
wavelength monitoring campaigns, such as the one reported 
in this paper, allow the detailed investigations that can reveal 
the processes leading to the formation of AGN broad-band 
spectra. 
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